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Abstract 

We calculate the Higgs decay rate of h — > Z^ by including the contributions from new scalars 
with arbitrary quantum numbers of the weak isospin (T) and hypercharge (Y) in the standard 
model. We find that our general formula for the decay rate of h — > Z7 matches with that for 
h — > 77 in the limit of mz = 0, but it is different from those in the literature. To illustrate our 
result, by taking the current 2a excess of the h — > 77 rate measured by the LHC, we examine the 
corresponding shift for the Z'j decay channel due to the new scalar. We show that the enhancement 
or reduction of the h — > Z7 rate only depends on the relative size of T and the absolute value of Y . 
Explicitly, we predict 0.76 < Rz^ = T(h — > Z7) /T 'sM(h — > Zj) < 2.05 by imposing the observed 
range of 1.5 < i? 77 = T(h — > Zjj)/TsM(h — > 77) < 2, which is independent of the number of 
multiplets and the couplings to the Higgs particle as long as the scalars are heavier than 200 GeV. 
This result provides a clear signature for the future LHC measurements to test physics beyond the 
standard model. 
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I. INTRODUCTION 



The discovery of the Higgs-like particle (h) by ATLAS [l| and CMS Q collaborations is 
a great triumph of particle physics. The following immediate question is whether this newly 
discovered particle is really the Higgs particle in the Standard Model (SM). To achieve this 
goal, one has to determine its couplings to all the SM particles. However, recent results from 
Higgs searches at the LHC have already shown a hint of new physics in the Higgs diphoton 
decay channel. The experimental data have shown that the diphoton decay rate is about 
1.5 — 2.0 times larger than the SM prediction, while the measurements in other channels 
including WW* and ZZ* agree with the SM Higgs properties. Due to the fact that the 
Higgs particle is electric neutral, its coupling to diphoton must be induced by some charged 
particles running in the loops, while in the SM, the W-boson and top-quark loops give the 
dominant contributions |3H7|. If the enhancement in the diphoton channel persists, it clearly 



indicates js-l26j that there must be some additional new charged particles which couple to 
the Higgs and mediate this decay process. Based on the electroweak (EW) SU(2)l x U(1)y 
gauge symmetry, these new particles should also contribute to the h — > Z-y channel and 



generica 
in Refs. 



ly lead to a shift in the Z'j decay width from the SM expectation. As pointed out 



27H33|. the simultaneous measurements of the 77 and Z7 channels at the LHC will 



provide us with valuable information about the structure of new physics. 

In this paper, we focus on the simple scenario in which extra contributions to the 77 and 
Z7 decay widths arise from some new scalars beyond the SM particle content. Some similar 



scenarios have already been explored in Refs. 



27 
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However, the formulae for the scalar 



27H32| are not consistent with 



contributions to the Z7 decay width used in the literature 
each other. Moreover, they cannot be reduced to the corresponding rates for h — > 77 by 
taking the limit of mz = and making the charge replacement. As a result, these numerical 
predictions for the Z7 decay width may not be reliable. Clearly, it is timely important to 
re-calculate the scalar contributions to the decay. 

This paper is organized as follows. In section [Til we show the detail calculations of the 
scalar contributions to the decay rate of h — > Z7. In section IHIl we give phenomenological 
analyses of the correlations between the 77 and Z7 channels in several specific sets as well 
as a generic one of the scalar multiplets. We summarize our results in section IIVI 
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FIG. 1. Feynman diagrams for a charged scalar particle (S) contributing to h — > Zj. 
II. CALCULATIONS OF SCALAR CONTRIBUTIONS TO h -> Z 7 

For a general scalar particle 5 with the third weak isospin charge T 3 and non-trivial 
electric charge Qs\^ the relevant Lagrangian involving S is given by 

C s = (D,S)\D»S) - m 2 s S^S - \ HS (rf H)(S* S) - \ HS (H^T a H)(S^T a S), (2.1) 

where H is the SM Higgs doublet; T a is the SU(2)l generators, and the covariant derivative 
involving 7 and Z gauge fields is defined as 

D^S = (0„ + ieQAp + iegzssZjS , (2.2) 

where g Z ss = (T^ — Qs s w)/( s wCw) with s w = sin6 w and c w = cos6 w (6 W being 
the Weinberg mixing angle). We note that the mass splittings at tree level among the 
components of the scalar S come from the last term in Eq. (12. ip . while they also receive 



loop-induced contributions of 0(100) MeV due to the gauge boson exchange diagrams [43] . 
which can be ignored in our discussion. On the other hand, the electroweak precision 
measurements, that is, the oblique parameters would constrain the mass differences 5m 



among the scalar multiplet [44| to be smaller than a few O(10) GeV [4l|, |45|. As a result, we 
assume the masses of the scalar multiplet are degenerate in this paper for simplicity. The 
trilinear coupling between Higgs and the scalar particle is —gtsshS^ S = —XnsvhS^ S after 
Higgs develops the vacuum expectation value, (H) = v/y/2. The amplitudes for the three 
Feynman diagrams shown in Fig. [1] can be easily written down: 



-iM {l) =e 2 Qg ZS s\ H sv^- d J 



d d l (21 + qY(2l + 2q + pfe z * (q)e** (p) 



(2Ti) d (I 2 - m 2 s ) [(I + q) 2 - m|] [(l + q + p) 2 - m|] ' 

. Ao\ 9^ ^ A-ri I a"' - qy \m - ^q - y, 

-% 



m 2 } ' 



■ KA (2) 2 n , 4-d f d " 1 (2l-qY(2l-2q-pyef l *(q)e^(p) 

/Jdj nv 



(2.3) 



1 We use the convention Q = T3 + Y/2 in this paper. 
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where we have taken the dimensional regularization. Note that by changing the integration 
variable I — > —I, —iAi^ can be transformed into the form that is exactly identical to 
— iM^\ so that we only need to the calculate 2(— iM^) + (— iM.^). With the standard 
procedures of the Feynman parametrization and the translation of the integral momentum 
/, we obtain: 

-i(2AfW + M {3) ) 

4 _ d f d d l {(21 + qy{2l + 2q + P y - [(I + q) 2 - m|]<T }ef (q)e** (p) 

e {2Q)gzss^HsVfji 



(2k) 



e 2 (2Q)g ZS sXHsv^- d T(3) I dxdy 



(I 2 - m|) [(/ + q) 2 - m|] [(Z + g + p) 2 - m|] 
(27r) f 



_(/2_ A 2)3 ^(/2_ A 2)2 



— 4:y(l — x — y)(p^q u — g^p • q) — (1 — 2x — 2y)(l — x — y)m z g 

(I 2 - A 2 ) 3 



C(?)C(p) 



(2.4) 



where A 2 = m| — x(l — x — y)m 2 z — y(l — x — y)va 2 h . In the second line of Eq. (12 .4p . we have 
used the following on-shell identities: 



p 2 = 0, q 2 



m 



(p + qf 



m 



p-e A *(p)=q-e z *(q)=0. 



(2.5) 



With the help of the dimensional regularization, it can be proved that in the second equation 
in Eq. (12 .4p . the first two terms exactly cancel with each other and only the last finite term 
is left, given by 

-i(2M {l) +M {3) ) 

4y(l — x — y)(p ll q u — g^ v p ■ q) 



ie 2 (2Q)g Z ss^Hsv 4 -d 
16tt 2 ^ 



J dxdy | 



[l-2x- 2y)(l - x - y)m\g 



in 



2 . AlV 



2 s -x(\-x- y)m 2 z -y(l-x- y)m\ 



(2.6) 



e?(?)ef(p). 



m| — x(l — x — y)m 2 z — y(l — x — y)m 2 l 
From this equation after the integration, the first term leads to 

ie 2 (2Q)g zss X HS v 



- i(2M {1) + M {3) ] 



16tt 2 



(pV - <Tp ■ q)e Z *(q)et*(p)A z ^(r s , X s ), (2.7) 



while the second one vanishes, where t$ = 4m s /m 2 L , Xs = 4m s /m 2 z and the loop function 
A Zl (x, y) is defined in Appendix [A] By combining the SM contributions and integrating the 
phase space of outgoing particles, we obtain the h — > Z7 decay width, 



T(h Zj) 



Q- 



512tt 3 h 



m 



in j 



2 \ 3 



A Zl - 



><HSV 



(2j2Qs-gzss)A z "(T S ,X< 

T 3 



(2.8) 
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with 



SM 



2 

V 



cot e w A^(r w , x w ) + 7Vc (2g t )(r 3 (t) -2Q^) Af K 



At) 



(2.9) 



where Tj = 4m-/m^, A, = 4m?/m| = W, £), and the summation is over the different isospin 
components in given SU(2)l multiplets. We point out that for the non-SM contribution to 
the decay rate of h — > Zj, there are extra factors of -1/2 in Refs. 29j,|31| and -1 in Ref. 32( in 
comparing with ours in Eq. (I2.8p . The different signs for the scalar contributions to h — > Z^y 
clearly lead to different results by taking limits to h — > 77. The modification in the partial 



decay width of h — > Z7 is then expressed in terms of the enhancement factor, given by 



R 



T{h -> Zj) 

TsmO -> zi) 



1 - NX HS - 



2^2Qs- 9zss 

T 3 



(2.10) 



where Tski{h — > Zj) represents the decay width in the SM and the factor N represents the 
degeneracy of the multiplet. For completeness and convenience for the later discussions, we 
also present the standard formula for the enhancement factor of the diphoton rate, 

2 



R 



77 



T(h -+ 77) 
FsM(h -> 77) 



l + iV^4 



T 3 



ATM + N c QtAJ] 2 (r t ) 



(2.11) 



where the loop functions Aj 7 (x) (j = 0, 1/2, 1) are defined in Appendix S It is straightfor- 
ward to show that our formula in Eq. ( 12 . 1 0[) for the Z7 decay can be retrieved to the one in 
Eq. (12. lip for the 77 mode when taking mz — > and making the replacement gzss Qs- If 
the mass splittings among components of the multiplet are taken into account, the deviations 
of the contributions in Eqs. (" 12 . 1 0[) and (12. lip are approximately proportional to Sm 2 /m 2 s . 
The modifications from the mass splittings are at percentage level and are negligible due to 
the electroweak precision measurements. Consequently, we will simply take Eqs. ( 12 . 1 0[) and 
(12. lip to analyze the qualitative behaviors in h — > 77 and h — > Zj decays. 



III. PHENOMENOLOGICAL ANALYSES OF CORRELATIONS BETWEEN THE 
77 AND Z7 CHANNELS 

With formulae in Eqs. (12 . 101) and (12. lip at hand, we are able to study the correlations 
between h — > 77 and h — > Z^ decay rates. According to the EW gauge symmetry, one 
expects that the particle responsible for the observed enhancement in the diphoton channel 

5 



could also lead to the modification of the Z7 rate. Furthermore, it is interesting to notice 
that the Z7 amplitude depends not only on the electric charges of particles running in the 
loops but also on their isospins, in contrast with the 77 case in which only electric charges 
can be probed. Consequently, a combined analysis of the decay widths in these two modes 
can provide us with important information on the EW charges of particles contributing to 
both processes. 



A. Singlet and Doublet Scalars 

In this subsection, we revisit two simple cases of singlet and doublet scalars under SU(2)l 
with the couplings: 

(I): J&=^*, (H): gf 83 ^ -±- (\-Q s ^) , (3.1) 

where Qs = 1, corresponding to (T, Y) = (0,2) and (1/2,1), respectively. The aim for 
this part of the investigation is twofold: (a) illustrating the features of the isospin and 



hypercharge dominations and (b) reexamining the results in Ref. [29|. In the two cases in 
Eq. (I3.ip . since only charged particles with Qs = 1 run in the loops, the enhancements 
for the diphoton rate are expected to be the same. However, due to the different isospin 
representations, those for the Z7 mode are distinct. From Fig. |2} we see that the Z7 rate in 
Case I (II) has a small suppression (enhancement) in comparison with the SM value. These 
features can be understood from the general formula in Eq. (I2.10p . In order to enhance the 
h — > 77 rate by a constructive interference, it is required that Xhs must be negative since 
^o 7 ( r s) anc I ^sm possess an opposite sign. In this case, the h — > Zj rate depends only on 
the sign of Qsdzss as ^hs < corresponds to both positive values of A^m and Ag 7 (rs, As) 
in the physically interesting mass regime. As a result, Rz-y is enhanced (suppressed) if 
QsQzss is positive (negative). Accordingly, it is easy to see that Qsgzss < (> 0) in Case 



I (II). Note that our result is opposite to that in Ref. [29|, in which the scalar contributions 
to h — > 77 and h — > Z7 are positive correlated in Case I and negative correlated in Case II 
due to the extra minus sign in the formula. 

The factor Qs ■ gzss can be examined in a general scalar multiplet. By using the identity 
Q = T 3 + Y/2, we have 

1 ( ~ Y\ /„ 2 Ys^y 



(a) (b) 
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FIG. 2. Contours of R~ n (solid lines) and Rz-y (dashed lines with the values in the yellow boxes) in 
rns-Xas plane for scalars of (a) singlet (T=0 and Y=2) and (b) doublet (T=l/2 and Y=l), where 
the solid (dashed) lines from top to bottom correspond to the values of (a) 1.25 (0.94), 1.50 (0.9), 
1.75 (0.86), 2.00 (0.82) and 2.25 (0.74); and (b) 1.25 (1.10), 1.50 (1.15), 1.75 (1.20), 2.00 (1.25) 
and 2.25 (1.35), respectively. 

When the isoweak charge T of the scalar multiplet is larger than its hypercharge Y, we find 
Qs ' 9zss > 0, which indicates the enhanced behavior in the Z7 channel compared with the 
SM prediction. In contrast, for the multiplet with 7 > T, we get Qs • gzss < 0, leading 
to a suppressed rate of h — > Z^y. Precisely, Cases I and II are the two simple but nontrivial 
archetypes of the two limits, respectively. 



B. Triplet Scalars with T = 1 and Y = 0, 2 



It is known that trip 
II seesaw model 
scalars, especia 



signatures 



53 



et scalars are introduced in many SM extensions, such as the type 
46H52|. Current LHC experiments also plan to search for these triplet 
ly the doubly charged particle in the Y = 2 case due to its interesting 



54J. Instead of a special model, we examine the generic feature of these 



triplets. In particular, we investigate their loop effects in the 77 and Z7 decay channels 
of the Higgs particle. For simplicity, we assume the multiplets of scalars possess the same 
masses and no mixings between the triplet scalars and the SM Higgs. For more detail 



analysis one can refer to Refs. 



38 



42l | . Our result is shown in Fig. [3j Clearly, the Z7 



rate is enhanced in both cases, which belong to the isospin domination class discussed in 
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FIG. 3. Contours of i? 77 (solid lines) and i?^ 7 (dashed lines) in ms-Xns plane for the triplet 
scalars of (a) (T,Y)= (1,0) and (b) (T,Y)=(1,2), where the solid (dashed) lines from top to bottom 
correspond to the values of (a) 1.25 (1.25), 1.50 (1.53), 1.75 (1.75), and 2.00 (2.00); and (b) 1.25 
(1.05), 1.50 (1.10), 1.75 (1.15), and 2.00 (1.20), respectively. 

the previous subsection. Furthermore, the enhancement in the Y = case is much larger 
than that in the Y = 2 one, which implies that increasing the value of Y will decrease the 
h — > Z"f decay rate for a given isospin multiplet. Moreover, in comparison with Case I of 
(T, Y) = (0,2) in Sec. III-A, T(h -»■ Z 7 ) for the case of (T, Y) = (1,2) for a given rate 
of h — > 77, say, 1.5 times the SM value, is slightly larger, showing that a large SU(2)l 
multiplet helps to increase the h — > Z7 rate. 



C. 5-plet Scalars with T = 2 and Y = 0, 2 

The weak SU(2)l 5-plet scalar particles are also of great phenomenological interest. For 
example, a recent study in Ref. 55| shows that with the 5-plet and Y = 2 scalar, the 
small Majorana neutrino masses can be generated at two-loop level. The 77 and Z 7 results 
are presented in Fig. HI Similar to the SU(2)l triplet cases, the 5-plet scalars give sizable 
enhancements in the h — > Z7 channel when the diphoton rate is lying in 1.5 to 2.0 times 
larger than the SM value. Explicitly, the enhancement factor Rz^ is also about 1.5 to 2.0 
(1.25 to 1.50) for the Y = (2) case. Again, the behaviors of these two cases can be expected 
by the rules of the isospin domination. We note that a bit less enhancement in the Y = 2 
case is due to the hypercharge cancellation. 
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FIG. 4. Contours of R~ n (solid lines) and Rz-y (dashed lines) in ms-Xns plane for the 5-plet 
scalars of (a) (T,Y)= (2,0) and (b) (T,Y)=(2,2), where the solid (dashed) lines from top to bottom 
correspond to the values of (a) 1.25 (1.25), 1.50 (1.53), 1.75 (1.75), and 2.00 (2.00); and (b) 1.25 
(1.10), 1.50 (1.15), 1.75 (1.25), and 2.00 (1.50), respectively. 

D. General Discussions 

In this subsection, we give a more comprehensive overview of how to use the combined 
analysis of the 77 and Z7 modes to investigate the nature of the new scalar multiplets. By 
assuming that the excess of h — > 77 is originated from a constructive interference between 
the SM and charged scalar loops, we obtain 

^ - 1 4[ArW) + N c Q*AT} 2 (n)] A z \r s , A*) Er 3 Qs • 9zss 



= 2.71 ^0 Vs, A 5 ) 4T(T + 1)^-3^4, 



(3-3) 



Al\r s ) 4T(T + 1) + 3F 2 
for a scalar with general charges of T and Y . From Eq. (13.31) . Rz^ increases with a higher 
value 2 of T since Aq 1 and A 1 ^ are both positive, whereas it decreases with a larger 
Furthermore, when T ^> \Y\, Rz^ is saturated to a maximal value with a fixed value of i? 77 , 
given by 

^/r7^t - 2m ^r^r' (3 - 4) 



The renormalization group running will drive the SU(2)l gauge coupling to violate perturbativity for 



2T+ 1 > 



The main point in this paper is to show the quantitative behaviors of the scalars with 



arbitrary EW quantum numbers. 
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FIG. 5. Rz^ as a function of the isospin T with i? 77 = 1.5 and ms = 300 GeV 
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FIG. 6. Rz~/ as a function of the hypercharge Y with i2 77 = 1.5 and mj = 300 GeV. 

which is definitely positive, implying that the Z'j decay width must be enhanced. This limit 
corresponds to the isospin domination discussed in Sec. III-A. In Fig. [5j we show 

Rz-y as a function of the isospin T with i? 77 = 1.5 and the scalar mass mj = 300 GeV. From 
the figure, we find that the saturated value of R™® x is 1.53. 

On the other hand, when we take the large hypercharge limit \Y\ ^> T, that is, the 
hypercharge domination case, we obtain the minimum value for the factor i?^ 7 . 



v^-i -°' 63 '^f' (3 ' 5) 

In Fig. El we show Rz-y as a function of the hapercharge Y with i? 77 = 1.5 and mj = 300 GeV. 
In the large-|F| limit, we get the asymptotic value of R^J 1 = 0.86. 

The above features become more transparent in the large scalar mass limit, m| — > oo. In 
this limit, we have that v4q 7 (ts, \ s ) = Aq 1 (t s )/2 = 1/6, and thus Eq. (13 .3p can be reduced 
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FIG. 7. Rz-y in the plane of general SU(2)l X U(V)y quantum numbers by fixing 72 77 = 1.5, where 
the solid and dashed lines represent the results of mj — > oo and mj = 200 GeV, respectively. 



to 



V^" 1 .^ 4T(T + l)c^-3FV 



V^-l" 1 - 36 ' 4T(T + 1) + 3V 2 • ^ 
In most of the physically interesting regime, Eq. (I3.6P is a very good approximation with 
the error of ten percent for the scalar mass above 200 GeV in which the ratio V^rg) * s 
numerically stable around 0.5. A remarkable feature shown in Eq. (I3.6P is that the h — > 
decay rate for the Higgs boson only depends on the EW quantum numbers assigned to 
the scalar multiplet for a given enhancement in the 77 channel. As an illustration, Fig. [7] 
shows the contours of Rz-y in terms of the general SU(2)l x U(1)y quantum numbers with 
i? 77 = 1.5 for both mj — > 00 and ms = 200 GeV, respectively. In the large mass regime, the 



2 



maximum and minimum values of Rz-y can be easily obtained, R^ x = [1 + 1.04(y / i? 77 — 1)] 
and i?™™ = [1— 0.31( A /'i? 77 — l)] 2 , corresponding to the isospin and hypercharge dominations, 
respectively. 

Finally, we can further demonstrate that the above bounds hold when we extend our 
discussion to arbitrary numbers of scalar multiplets with different quantum numbers, masses 
rrii, Higgs couplings A i5 and the degeneracies Aj. The formula in Eq. (13. 3 p can be generalized 
to: 

_ 136 gj j^jkgj + + Iffi - 3V, 2 ^] 

v^- 1 ' ' E,^S( 2 ^ + 1 )[ 4T ^ T ^ + 1 ) + 3 ^ 2 ] 

where we have taken > 200 GeV. Note that Eq. (13.71) is valid only for the new scalar 

11 



contributions interfering with the SM part constructively. This assumption follows by two 
indications: (i) Aj should be negative in sign, so that we can replace A; with its absolute 
value in Eq. ( 13. 7ft ; and (ii) contributions to the amplitudes of the h — > 77 and Z7 modes 
must be smaller than the corresponding SM parts in size to eliminate the sign ambiguity 
when we take the square roots of _R 77 and i?z T We now try to derive a generic bound for 
Rz-y- From Eq. (13. 7\i and inequalities 



- 0.31 [4Tj(Tj + 1) + 3Y 2 ] < 4Ti(Ti + l)c 2 w - ZY?s 2 w < 1.04 [4T,(T, + 1) + 31?] (3.8) 



From the general bound in Eq. (13.91) . it is clear that one can deduce the allowed range of 
the Z7 decay width for a given value of _R 77 . Since the current ATLAS and CMS data have 
indicated the range 1.5 < R~ n < 2.0 for the diphoton mode, one predicts that 



for the Z7 mode, where both the upper and lower limits arise from i? 77 = 2.0. We emphasize 
that the result in Eq. (13 . 1 0[) does not rely on the details of the extra scalars, such as the 
quantum numbers and their couplings to the Higgs particle. Thus, if the future-measured 
diphoton rate is still enhanced about 1.5 — 2.0 times over the SM value, the Z7 decay rate 
should be in the range of 0.76 to 2.05. Otherwise, our simple scenario discussed in this paper 
will be excluded. 

IV. SUMMARY 

We have studied the decay width of the h — > mode with extra scalar multiplets of 
arbitrary weak isospins and hypercharges. Due to the SM gauge symmetry, these new scalars 
which result in the excess of the h — > 77 decay rate would generically shift the h — > 
decay rate as well. The combined analysis of these two modes can provide us with valuable 
information about the new physics structure. Regarding the ambiguity of the formulae for 
Y(h — > Z'-f) in the literature, we have revisited some simple extensions to the SM in which 
new scalar multiplets are introduced. We have found that T(h — > Z7) only depends on the 



for a single multiplet, we find 




(3.9) 



0.76 < R z ^ < 2.05 



(3.10) 
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relative size of the isospin T and the absolute value of the hypercharge Y. In particular, we 
have shown that the enhancement factor i?^ 7 is a monotonically increasing function of T and 
a monotonically decreasing one of |V| with the fixed value of i? 77 . This observation enables 
us to predict that 0.76 < Rz~, < 2.05 by imposing the observed range of 1.5 < i? 77 < 2 if 
the scalars are heavier than 200 GeV. Note that this range is irrelevant to the number of the 
scalar multiplets, their representations, and the couplings to the Higgs particle. Our results 
on the h —> Zj decay clearly can be tested at the LHC. 



Appendix A: Definition of Loop Functions 

The functions related to h — > 77 and h —> Z7 are defined as follows: 

AP(x) = -x 2 [2x~ 2 + 3X" 1 + 3(2x" 1 - (Ala) 

AJ] 2 (x) = 2x 2 [x~ 1 + (or 1 - l)f(x-% (Alb) 

Al\x) = -x 2 [x~ l -f{x~% (Ale) 

Af 7 (x, y) = 4(3 - tan 2 6 w )I 2 {x, y) + [(1 + 2X' 1 ) tan 2 9 W - (5 + 2x~ 1 )]/i(a;, y), (Aid) 

Ay 2 (x, y) = h(x, y) - h(x, y), (Ale) 

Ap(x,y) = I 1 {x,y), (Alf) 

where 

= K^y) + W^ 1 ^ ~ /(!rl)1 + T^V? [g[X ' l) ~ g(rl% 

h(x,y) = -^r^lf^ 1 ) - /(»"')]■ (A2b) 

For a Higgs mass smaller than twice of that of the loop particle, i.e. rrih < 2mj 00 p, we have 

f(x) = arcsin 2 y/x, (A3a) 



g(x) = Vx 1 — 1 arcsin yfx. (A3b) 



ACKNOWLEDGMENTS 



The work was supported in part by National Center for Theoretical Science, National Sci- 
ence Council (NSC-98-2112-M-007-008-MY3 and NSC-101-2112-M-007-006-MY3) and Na- 

13 



tional Tsing-Hua University (102N1087E1), Taiwan, R.O.C. 



[1] G. Aad et at. [ATLAS Collaboration], Phys. Lett. B 716, 1 (2012) [arXiv:1207.7214 [hep-ex]]. 
[2] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716, 30 (2012) [arXiv: 1207.7235 
[hep-ex]]. 

[3] M. A. Shifman, A. I. Vainshtein, M. B. Voloshin, and V.I. Zakharov, Sov. J. Nucl. Phys. 30, 

711 (1979) [Yad. Fiz. 30, 1368 (1979)]. 
[4] J. R. Ellis, M. K. Gaillard, and D. V. Nanopoulos, Nucl. Phys. B 106 (1976) 292. 
[5] B. L. Ioffe and V. A. Khoze, Sov. J. Part. Nucl. 9, 50 (1978) [Fiz. Elem. Chast. Atom. Yadra 

9, 118 (1978)]. 

[6] T. G. Rizzo, Phys. Rev. D 22, 178 (1980) [Addendum-ibid. D 22, 1824 (1980)]. 
[7] L. Resnick, M. K. Sundaresan and P. J. S. Watson, Phys. Rev. D 8, 172 (1973). 
[8] F. J. Petriello, JHEP 0205, 003 (2002). 

[9] T. Han, H. E. Logan, B. McElrath and L. -T. Wang, Phys. Lett. B 563, 191 (2003) [Erratum- 

ibid. B 603, 257 (2004)] [hep-ph/0302188]. 
[10] C. -R. Chen, K. Tobe and C. -P. Yuan, Phys. Lett. B 640, 263 (2006) [hep-ph/0602211]. 
[11] G. Cacciapaglia, A. Deandrea and J. Llodra-Perez, JHEP 0906, 054 (2009) [arXiv:0901.0927 

[hep-ph]]. 

[12] J. Cao, Z. Heng, T. Liu and J. M. Yang, Phys. Lett. B 703, 462 (2011) [arXiv: 1103.0631 
[hep-ph]]. 

[13] S. Dawson and E. Furlan, Phys. Rev. D 86, 015021 (2012) [arXiv: 1205.4733 [hep-ph]]. 
[14] T. Kitahara, JHEP 1211, 021 (2012) [arXiv: 1208.4792 [hep-ph]]. 

[15] A. Joglekar, P. Schwaller and C. E. M. Wagner, JHEP 1212, 064 (2012) [arXiv: 1207.4235 
[hep-ph]]; 

[16] M. Hashimoto and V. A. Miransky, Phys. Rev. D 86, 095018 (2012) [arXiv: 1208. 1305 [hep- 
ph]]. 

[17] A. Delgado, G. Nardini and M. Quiros, Phys. Rev. D 86, 115010 (2012) [arXiv: 1207.6596 
[hep-ph]]. 

[18] H. An, T. Liu and L. -T. Wang, Phys. Rev. D 86, 075030 (2012) [arXiv:1207.2473 [hep-ph]]. 

11 



[19] N. Arkani-Hamed, K. Blum, R. T. D'Agnolo and J. Fan, JHEP 1301, 149 (2013) 

[arXiv:1207.4482 [hep-ph]]. 
[20] L. G. Almeida, E. Bertuzzo, P. A. N. Machado and R. Z. Funchal, JHEP 1211, 085 (2012) 

[arXiv: 1207.5254 [hep-ph]]. 
[21] E. Bertuzzo, P. A. N. Machado and R. Zukanovich Funchal, JHEP 1302, 086 (2013) 

[arXiv:1209.6359 [hep-ph]]. 
[22] G. Moreau, Phys. Rev. D 87, 015027 (2013) [arXiv: 1210.3977 [hep-ph]]. 
[23] M. Chala, JHEP 1301, 122 (2013) [arXiv: 1210.6208 [hep-ph]]. 

[24] S. Dawson, E. Furlan and I. Lewis, Phys. Rev. D 87, 014007 (2013) [arXiv:1210.6663 [hep-ph]]. 
[25] K. Choi, S. H. Im, K. S. Jeong and M. Yamaguchi, arXiv:1211.0875 [hep-ph]. 
[26] R. Huo, G. Lee, A. M. Thalapillil and C. E. M. Wagner, arXiv:1212.0560 [hep-ph]. 
[27] A. Djouadi, V. Driesen, W. Hollik and A. Kraft, Eur. Phys. J. C 1, 163 (1998) [hep- 
ph/9701342]. 

[28] A. Djouadi, Phys. Rept. 457, 1 (2008); B. A. Kniehl, Phys. Rept. 240, 211 (1994). 
[29] M. Carena, I. Low and C. E. M. Wagner, JHEP 1208, 060 (2012); 
[30] C. W. Chiang and K. Yagyu, arXiv:1207.1065 [hep-ph]. 
[31] I. Picek and B. Radovcic, arXiv: 1210.6449 [hep-ph]. 
[32] W. C. Huang and A. Urbano, arXiv:1212.1399 [hep-ph]. 

[33] G. n. Li, G. Guo, B. Ren, Y. J. Zheng and X. G. He, arXiv: 1212.5528 [hep-ph]. 

[34] C. Han, N. Liu, L. Wu, J. M. Yang and Y. Zhang, arXiv:1212.6728 [hep-ph]. 

[35] B. Batell, S. Gori and L. T. Wang, JHEP 1206, 172 (2012) [arXiv:1112.5180 [hep-ph]]; 

[36] S. Baek, P. Ko and E. Senaha, arXiv:1209.1685 [hep-ph]. 

[37] W. F. Chang, J. N. Ng and J. M. S. Wu, Phys. Rev. D 86, 033003 (2012) [arXiv: 1206.5047 
[hep-ph]]. 

[38] A. Melfo, M. Nemevsek, F. Nesti, G. Senjanovic and Y. Zhang, Phys. Rev. D 85, 055018 

(2012); [arXiv:1108.4416 [hep-ph]]. 
[39] A. Arhrib, R. Benbrik, M. Chabab, G. Moultaka and L. Rahili, JHEP 1204, 136 (2012); 

[arXiv:1112.5453 [hep-ph]]. 
[40] A. G. Akeroyd and S. Moretti, Phys. Rev. D 86, 035015 (2012); [arXiv: 1206.0535 [hep-ph]]. 
[41] E. J. Chun, H. M. Lee and P. Sharma, JHEP 1211, 106 (2012); [arXiv: 1209. 1303 [hep-ph]]. 
[42] P. S. B. Dev, D. K. Ghosh, N. Okada and I. Saha, arXiv:1301.3453 [hep-ph]. 

15 



[43] M. Cirelli, N. Fornengo and A. Strumia, Nucl. Phys. B 753, 178 (2006). 

[44] L. Lavoura and L. -F. Li, Phys. Rev. D 49, 1409 (1994) [hep-ph/9309262]. 

[45] M. Baak, M. Goebel, J. Haller, A. Hoecker, D. Kennedy, R. Kogler, K. Moenig and M. Schott 

et al, Eur. Phys. J. C 72, 2205 (2012) [arXiv:1209.2716 [hep-ph]]. 
[46] M. Magg and C. Wetterich, Phys. Lett. B 94, 61 (1980). 
[47] J. Schechter and J. W. F. Valle, Phys. Rev. D 22, 2227 (1980). 
[48] T. P. Cheng and L. -F. Li, Phys. Rev. D 22, 2860 (1980). 
[49] G. B. Gelmini and M. Roncadelli, Phys. Lett. B 99, 411 (1981). 
[50] G. Lazarides, Q. Shafi and C. Wetterich, Nucl. Phys. B 181, 287 (1981). 
[51] R. N. Mohapatra and G. Senjanovic, Phys. Rev. D 23, 165 (1981). 
[52] J. Schechter and J. W. F. Valle, Phys. Rev. D 25, 774 (1982). 

[53] S. Chatrchyan et al. [CMS Collaboration], Eur. Phys. J. C 72, 2189 (2012) [arXiv: 1207.2666 
[hep-ex]]. 

[54] G. Aad et al. [ATLAS Collaboration], arXiv: 1210.5070 [hep-ex]. 

[55] C. S. Chen, C. Q. Geng, D. Huang and L. H. Tsai, arXiv: 12 12.6208 [hep-ph]. 



16 



